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ABSTRACT 


ANALYSIS OF INELASTIC ALPHA PARTICLE SCATTERING 
USING DIRECT INTERACTION THEORY 


by 


Herbert Frank Swanson Jre 


The theory of the direct interaction is discussed 
and its applicability to inelastic alpha particle 
scattering to low lying states is shown. Two results 
of this theory are developed; (1) the diffraction model 
of Blair and (2) the distorted wave Born approximation, 
where the phenomenological approximation to the optical 
wave function suggested. by McCarthy and Pursey is used. 
Fits to Ароо( оќ, 6) Аро data are shown for 32.8 and 
41.0 Mev(Lab) alpha particles. In the diffraction model 
analysis, К. = 6.45 Fermis for both energies. | 
ARS = „82 Fermis for the lower energy data and .71 
Permis for the higher energy data. ZAsRo = .67 and „56 
Fermis for the lower and higher energy, respectively. 
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CHAPTER I 


THE DIRECT INTERACTION 


There are two basic mechanisms by which a nuclear 
reaction can proceed; (1) the direct interaction, and 
(2) the formation of a compound nucleus. The theoretical 
angular distribution predicted by: each mechanism is 
different. To see the nature of these two processes, 
consider the incident particle to be a nucleon. Thus 
the independent-particle shell model is a good 
approximation to the nuclear structure. Here the many 
body problem is replaced by an effective single particle 
potential which describes the force acting on each 
particle, this potential being different for different 
shell model configurations. 

As a particle enters the nucleus, it experiences 
a potential described by the other A nucleons in their 
ground state.configuration. The projectile may then 
simply enter the nucleus, be deflected by the nuclear 
potential, and emerge again at a different angle but 
with the same energy in the center of mass system. 
This is called direct elastic scattering. If a 
reaction ensues, the energy transferred to the nucleus 
is determined by the interaction of the incoming 


nucleon with one in some shell configuration and by 





the nuclear density. This later fact modifies the 
problem as neither of the colliding particles can have 
an energy below that of the ground state energy of the 
nucleus. This is because the hole left by the struck 
nucleon is the only level available below the ground 
state energy. The case where the incident particle 
replaces the struck nucleon is called ος scattering 
and will not be considered here. Those nucleons excited 
to high level states should be described by a different 
potential than those excited to low lying states, 

hence the effective potential is a function of the 
bombarding energy. If either of the colliding nucleons 
has an energy greater [πε its separation energy, it 
may leave the nucleus without further reaction save 
deflection by the nuclear potential. This is described 
as a direct interaction and usually occurs in a time 
comparable with the time it takes the projectile to 
traverse the nucleus. If, on the other hand, neither 
of the nuclei has energy sufficient for separation, 
doth will undergo further collisions eventially 
spreading their energy over the whole nucleus. For a 
time, the nuclear state will grow increasingly complex 
until it reaches a statistical equilibrium. This state 
will be a complicated admixture of shell-model and 

. collective configurations involving excitation of the 


many degrees of freedom of the nucleus. Eventually 





enough energy will be concentrated in a particle or 
group of particles to allow escape from the nucleus 
and upon emission the nucleus will decay into some 
new excited state. This process is referred to as 
scattering via the formation of a compound nucleus 
and is usually characterized by the long time necessary 
for its completion, which is about a thousand times 
that for the direct interaction. Since these two 
processes exhibit widely different cross sections, a 
method for Separation into the direct interaction 
and coitum nucleus components is desirable. No 
suitable technique, however, has been found. Instead 
reactions are chosen which are predominently either 
direct interaction or compound nucleus types and where 
errors made in neglecting the other are small. 

Direct interaction processes involve only a few 
of the mere decrees of freedom of the nucleus. The 
minimum number, in fact, would be just those necessary 
to specify the ingoing and outgoing channels of the 
reaction. The concept of the direct interaction can be 
extended to any projectile if it excites only one degree 
of freedom in the nucleus specifically a collective 
degree of freedom. For example, it would certainly be 
& direct interaction if the projectile bounced off a 
. Spheroidal nucleus, leaving it rotating. These 


rotational states are typically low lying states. 





Compound nucleus reactions, however, involve 
many degrees of freedom due to the complicated nature 
of the nuclear excited state. It would then be 
fortunate to have a detector capable of indicating 
the number of degrees of freedom in the wave functions 
of the intercepted particles, but alas, none has ever 
been found. In a typical compound nucleus reaction 
the probability of a particle being emitted leaving 
the nucleus in a low lying state is greater than those 
leaving the nucleus in the higher excited states. 
This is due to the difficulty of barrier penetration 
which inhibits the emission of particles with energies 
just above the separation energy. Figure (1-1) shows 
a rough plot of the number of particles emitted vs. 
the excitation energy of the nucleus, where ες is the 


ground state energy. 
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The vertical lines show the discreet spectrum. The 


detector, however, does not possess perfect resolution 


E 





and the energy spectrum it sees is averaged over an 
interval AE. The solid curve shows such an averaging. 
Thus it is seen that although the compound nucleus 
would rather emit particles leaving tne energy of tne 
nucleus near its ground state, the majority of the flux 
occurs in the higher lying states. 

Hence it is seen that inelastic scattering to 
low lying states should go predominantly by the direct 
interaction process. This theory will be developed 


in the following chapters. 





CHAPTER 11 
THE SCATTERING CROSS SECTION 


Introduction 


A useful concept in relating theory to experiment 
is the differential cross section. In the usual picture 
of scattering, a collimated beam of particles impinges 
upon a target. The number of particles per second, an, 
scattered into an element of solid angle is proportional 


to the incident flux N and the element of solid angle dí. 


dN = c-(eq) N ds er) 
The constant of proportionality 0(6,) has dimensions 

of area and is called the differential cross section. 

It can be a function of O and 0 the angular coordinates 

of AN, but is usually only a function of O taking the 
polar axis as the beam direction. The cross section is 
related to Schrodinger's equation in the following manner. 
At distances far from the scattering center, the wave 
function describing the scattered particles can be 
written as the sum of a plane wave corresponding to the 


incident beam and an out going spherical wave due to the 





scattering. 
С. 
Мае? а у, 
(2-2) 
Y CLE 
OCAT = L(09) 8 
n^ 


The function Íte 6) is known as the scattering anplitude 
and contains all of the angular dependence of the’ 
scattered beam. | 

The quantum mechanical expression for the scattered 


ii 
flux is given by 


ames et sns УХ Scar — us] Е 


where me is the reduced mass of the scattered system. 
Applying this to (2-2), an expression is obtained 
giving the number of particles being scattered / unit 


area / sec. 


ο 


Sar me pz (2-4) 
where, ín the asymptotic region considered, only the 
radial component contributes significantly to the flux. 

The number of sc se s particles passing through 
an element of area GA of a sphere of radius r located 


concentric with the scattering center is just dN. 


~ 





dN = Cser) dA = ES (2-5) 
where do. = dA 
pe 


The incident flux is computed using an expression 


ПЕЕ 
, Similar to (2-3), where the incident wave is € 


Sina = nes S (2-6) 


I 
The quantity dN is then divided by the incident flux 
and the resulting expression compared with (2-1) to 


Obtain the scattering cross section; 


seos (B () [Rea em 


Im the above intuitive argument, it should be 
noted that any interference between the incident and 
scattered waves has been ignored. It would thus appear 
that (2-7) is inaccurate and of limited usefulness. 

The incident beam however, has a macroscopically small 
cross sectional area and thus would only interfere with 
the scattered wave at both very small and very large 
scattering angles. Since the density of particles in 
these beams is small, they will not greatly influence 





One another even in these cases. For amore rigorous 
treatment, the reader is directed to reference / 2 7, 
where it will be noted that the result is the same, 

It is the purpose of the remaining sections of 
this chapter to obtain expressions ror the scattoring 
amplitude, using the Schrodinger equation with an 
appropriate potential. As this is a study of the 
direct interaction, the Hamiltonian will exclude all 
flux that would go into populating the states of the ` 
compound nucleus and would itself describe only those 
few internal degrees of freedom pertinent to the direct 
interaction. Typically, it would describe the motions 
of the one or two nucleons mich momentum had been 
transferred. 

| 


The Schrodinger Equation 


Consider the Shrodinger equation in the center 
of mass system. The subscript f means "with respect 


to the final channel f ." 


Er Ut) VIE) $ Z w —_ СА (2-8) 


The first two terms on the left are the kinetic energy 
and the optical interaction potential of the separating 
fragments, and AE) is the Hamiltonian describing 


their internal motions. VCH 5) is the residual 





LO 


interaction energy, with respect to the final state 
channel, which gives rise to transitions to excited 
states of the residual nucleus. The optical potential 


is of the form 


2 e° ντ 
Е zs (+ e ^7 9/4) 


where the term on the left is tho Coulomb potential 

for a point charge Z and an incident particle of 

charge 4 . The term on the rignt is the phenomenological 
complex nuclear potential first introduced by Saxon 

and thus bears his name. VWhile it will suffice here 

to give a plausibility argument for the complex nature 
of this potential, the reader is referred to / 3 / and 
Га У for a theoretical justification. The study of 
dispersion theory in Davee optics has shown that 
absorption of light in a cloudy or semi-transparent 
material can be represented mathematically by the 
imaginary term of a complex index of refraction. In 
the study of nuclear reactions, the phase of a particle 
wave impingent on a potential is shifted as it passes 
through the region of the potential and it is known 
that the scattering amplitude is completely describable 
in terms of these phase shifts.? Completely analogous 
to the optical case, 1f these phase shifts are complex, 


the imaginary part would produce absorption of part of 





ae 


the incident beam. It is then the imaginary part 

of the complex potential which gives rise to absorption. 
This term literally absorbs out all flux which is not 
described by the potential VG), such as the compound 
nucleus flux, thus eliminating it from the scattered 
beam. The denominator gives the radial dependence of 
the potential combining its short range nature with a 
diffuse edge. It will be mentioned in passing that 
many variations on this appear in the literature. 
However, they often contain so many parameters that the 
question arises as to whether or not the resultant fits 
contain any physics or are just attempts at I 
Piecing. The variables zc (or.&. in the incident channel) 
are those few internal degrees of freedom considered. 
The Hamiltonian K ( С.) is assumed to have eigen 


functions Vel Se) which are solutions of tne equation 


HE) VE) = hve) (2-9) 


To sOlve equation (2-8), it is first written 
with the two potential terms isolated. 


(Ev zy sue ea 


eme 


To remove the 5 dependence, the equation is multiplied 
by the complex conjugate of Vc Cz y and integrated 


over the internal variables. 





12 


- А-у? = | weds + fi FRE) IE 


Ze 
(2-11) 
% 
=- Jy (αν) 48, 
Using (2-9) and the Hermitean nature of HE, 
equation (2-11) can be written, using Dirac's 
"bra-ket" notation for the integrals as 
2, [2 LM 
(y +), = dx (ve |u*v]^?» (2-12) 
oM. = (Ve Cg) | Ὃν | 8 
E re = 


We is the coefficient in an expansion of ^l^ into a 
set of eigen functions vc and can be thought of 
as a projection of Y on the basis vector for the final 


The homogeneous part of Equation (2-12) 
6 


enannel Ur", 
is the Helmholtz wave equation with solutions 


ф = Bode С 


General Integral Equation Solution 


Equation (2-12) is transformed into an integral 
equation through the use of a Green's function. It 


can be shown that 


V 


ο... 


is a solution to (2-12), where the Green's function 


Ks): is a solution to the equation 


ν΄: Е.) Kite, £2) = Á (r.- z) (2-16) 


as long as Ppk re $ (2-16) is just the wave equation 


mentioned earlier with outgoing solution 


° = al ke Се 
This is expanded in a series of Legendre polynomials” 
. co и B Y 
eke = 2 0) зе) Т 
20 


Using the Legendre addition theorem", (2-17) can be 


expressed in the space fixed coordinate system. 





РА = “Ἐν („зү өө (2-18) 


{»о мя-Ё 


The primed and unprimed variables are the angular 
coordinates of kp and rp respectively. At large 
` distances from the origin, the asymptotic form 

of Jg (kr) may be used. 


4) (К.Е) Ж Sin Сага "x (2-19) 
NE с СМ", DEN А: 


0-2 ima-2 


We proceed by developing a series solution for the 
Green's function. All "f" subscripts will be dropped 
until they again become necessary. The three 
dimensional Dirac delta function expressed in spherical 


polar coordinates is? 


Ыс-с) = - 1. 90-0) 5(0:0-0:0) 500-0) (2-21) 


The spherical harmonics form a complete set over the 


surface of a sphere and a mathematical statement of 


this is the following "closure relation"? 





to 


ó(Cse- αμ dH (oleo (2-22) 


AD Mz, 


The total delta function can then be expressed in'terms 
of this series. The closure relation suggests the 


following expansion for the Green's function. 


Deo wa = 


(ος!) = б. А (οφ) аон) Yale® (2-23) 


This is now inserted into the defining equation, 


(2-16), where the Laplacian in spherical coordinates is 





.4 909, 1 C 2% a 
VE u$ + 2 i pioa 59: (2-20) 


Mo (2-16) becomes 


E. 5 A EE: (60) y Pa ΠῚ 


mal 
* ds [oe oro 2r) + с ahe Me 
-- FEO es 


The spherical harmonics are a solution to the following 


differential equation. 





| 96 
Sino o6 











"NC 
sirva e + 40+) ү, = (у 


This is solved for the second set of square brackets 


and the two series are equated term by term. 


A (eq) E Y, (ed) (2-26) 
NONO) xs (E “(εὐ ra rr) = Köln) (2-27) 


E 91079) 2 raten) (2-28) 


The radial Green's function is a solution to 


29694 (2 slit) Geria- 4 ölr-r) (2-29) 


Consider the homogeneous part of this equation. 


а Fler) 4 (E -4U) £p) SO. . (2-30) 
pe tok A | 


As this is a second order differential equation, 
it has two independent solutions. The radial Green's 
function, ©, P), is built from these two solutions 


by the folloving prescription. 





27 


Region I Recon τα 





FIGURE (2-1) 


Figure (2-1) shows a plot of some general G( sr) 
as a function of r. It must be a solution to (2-30) 
in both regions, satisfying the boundary condition 
at the origin in region I and the boundary condition 
at IntCinity in ropglon II, wile it its required that 
the magnitude of T) (sr) be continuous across the 
boundary at r=r’, there will be a discontinuity in 
the slope. The magnitude of this discontinuity can be 
obtained by integrating equation (2-29) from r-e tor'+e 

and taking the limit as €-vO. 


reé | че 
Limit / d ($667) + / EZ 4 (2) ) гдан 
eT Re τ ia 
i rE- 
== |2 SY di 
= 
ra 


d N al, у El. = 


ptc 
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The second integral contributes nothing, as Gılr,r’) 
is continuous across the discontinuity at r'e Two 
solutions to equation (2-30) will be obtained. 


Introduce the substitution 


f (kr) = /kr (ЕК) 
This is inserted into the homogeneous equation (2-50), 
yielding 


фе. | Д7 (а LE 
82 + -L 92 +(Е чы) = 0 (2-32) 


pen 


This is Bessel's equation of index ¿el and 


Z (kr) = κ 
Kis 
Nultiply through by the constant /Z 


$9 = kr E. J Ско) = με ύέᾳ (2-33) 


where £n is the spherical Bessel function. 
Consideration of equation (2-17) shows this to be 

kr times the radial wave function. Here, 44 (en is 
regular at the origin and thus satisfies the boundary 
condition in region I. The following relation exists 


between the spherical Bessel and Hankel functions.++ 


Ay ber) = A (hy len + (ЕВ) ^— (2-34). 
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() (2) 
The functions h, and y are, respectively, tne 


spherical Hankel functions of type one and type two. 
They are both solutions to equation (2-30). 

азупрьб ол и as an outgoing spherical wave 
and id an incoming ono. Аз г approaches infinity, 
there can be only outgoing scattered vaves. Hence 

a solution satisfying the boundary condition in 

τεσ!" ΤΙ Auer). The Green!s function is writven 


C, kr tier] rar! : 
G =) = Ὃν (2-35) 
C. kr? у бек) Pere | 


subject to the condition imposed by (2-31). 


d [e] a qero ]= as 


har ra 
The Wronskian relation for these two functions can be 


used to evaluate the constants Cj and Co. 
Sturm-Liouville theoryl^ gives the Wronskian of two 


independent solutions to equation (2-30) as 
| q | 
W (y; ter) ЕЮ) = A . (2-86) 


where A is a constant. 





20 


Wehe) = (634 i, 6f) 


ET A S ee 44) | (2-37) 


As this is true for ail Ὁ, the asymptotic forms 
of these expressions may be used in computation. 
The asymptotic form of per) vas already given 
in (2-19) and for large r 


, 2 
(1) | CC e - =) 
eel (2-38) 
hy Cer) — n 
Calculation of' the Wronskian yields 
А = К 
A RIOS of (2-31a) with (2-37) shows 
(1) 
C, = My (ler) 
eam t 
б. = фе) 
The Green's function becomes 
Gu, ` А 
Соб") = кеаби) h, (kh) C 


Y is the smaller of r and r’ 


F. is the greater of r and г” 





ll 


Combining equations (2-23), (2-26), (2-28), and 
(2-39), the total p Me function is obtained. 


Kg) = EY» рва ес 


(so ma 
x AG &) Yes 


Equation (2-23a) is then used in (2-15) to give 


(2-232) 


the solution for ^p 


The General Scattering Amplitude 


As was discussed in the introduction of this 
chapter, the scattering amplitude is the coefficient 
of the outgoing scattered spherical wave. Equation 
(2-15) will be used to find the asymptotic form 
of jr. Equation (2-23a) can be simplified for large р, 
г being 7» r; where " can be replaced by its 
asymptotic form given in (2-39). The Green's function 


becomes 


He 


Κίς, τ) Ξ 5 co^, (2, AY Ce fed (2-230) 


ΟΝ Ma 


l 
where the primed angular coordinates refer to Г. апа 
the unprimed ones to which is taken to be along 
the vector poen Consider the expression in (2-18). 


Taking the complex conjugate yields 
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e Es пер 4 52i: cc) “2 E Ced 


Azo Mak 


where the primed coordinates refer to i and the 
unprimed ones to e - | 


es = Ὁ 


Mods.) = ine 


Using equation (2- 15) 


ec ee 


yore Let Πα. | {6 “γαενγεδ-ᾶς 
A 


The plane wave e ““is a solution to the. homogeneous 
part of (2-12) and can always be added to the scattered 
solution. Thé scattering amplitude is then equated 


Z ο 
with the coefficient of = e ee 
РА 


Ко - = THe С Ls y | (2-40) 


where the Dirac notation is again used to indicate 


integration over all internal and external variables. 


The Elastic Scattering Amplitude . 


In a description of only the elastic component 
of the scattering, the potential We E) which 


describes transitions to excited states is set equal 





2Š 


to zero and the interpretation of the imaginary part 
of the complex potential is broadened to include 
absorption of the entire inelastic flux. The initial 
and final channels are then the same and the outgoing 


wave function Wis 
| ¢ 
Ф = vw, Yi, и (2-41) 


where νε( Σχ) and Y Ke, 2 are, respectively, the 
internal and elastic scattering wave functions for 


the initial channel. satisfies the equation 
H 2 | J Y 
CHV T и) (kg ti) = EX Kg tz) 


The elastic scattering amplitude is obtained from 
(2-40) 


fled = es és | ake > (c 


Inelastic Scattering - Distorted Waves 


The direct interaction inelastic scattering 
component is to be isolated from the scattered wave 
and. using νο ο ος in the first part of 
this chapter, the inelastic scattering amplitude will 
be found. Equation (2-8) is re-written with the 


potential which describes the inelastic transition d 


= 





to 
e 


isolated. 


Kyu. ХЕ) -= Y = VU (E) Y 


vene 


Using (2-9) and the Hermitean nature or ( =.), 
the E dependence is removed in the same manner as 


for (2-12) 


(yy + 4 )^p. = BEE | VER, EJ Y) (2-45) 


Us ' SZ C Cr) | (2-44) 


, | 
Ko and ^P. аге discussed after equation (2-12) and will 


be reiterated here for the readers! convenience. 


V = «Чо [у | 
К. = E las h) 


The homogeneous part of (2-43) has as one of its 
eigen solutions, the elastic scattering wave function 
discussed in relation to equation (2-41). A series 
solution corresponding to an outgoing wave will be 
obtained. Consider the homogeneous equation with 


the f subscript dropped as before. 


-U+ к) © (2-45) 
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Tne form of the desired solution is such as to reduce 
ES expression (2-18) when the optical potential is 
removed, as certainly (2-45) reduces to the Helmholtz 
wave equation. As the optical potential is a function 
of only the radial variable, Xen will have 
azimuthal symmetry if k is assumed to be along tne 
z-axis. (Note expression (2-17)). The discussion 
following Equation (2-8) gave the optical potential 
to be of the form | 

N anod EF K a 
A plausibility argument for the complex form of the 
nuclear potential was given as well as is expression, 
typical of those appearing in the literature. No 
explicit form of Vy is necessary in the following 
calculations. However, it will be assumed that the 
nuclear forces are short ranged and thus Vy will 
rapidly approach zero for distances greater than the 
nuclear radius. Equation (2-45) becomes 


[v^ «E - £85 — ere as po Ó . (2-452) 


zma 2. | 
where = οἱ - E = Kr (2-46) 
| vu” 
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Coulomb Wave Functions 


Consider the special case of zero nuclear 
potential. This situation would be encountered if 
the impingent beam lacked sufficient energy to 
Overcome the coulomb barrier and penetrate into the 
nuclear material. The eigen solutions become coulomb 


wave functions, satisfying the equation 
ΚΩ Е - «η Q Ce) = ` (2-47) 


\ 
The second chapter of reference / 3 / is an excellent 
source for information regarding solutions of tnis 
equation. It turns out that (2-47) has a regular 


solution of the form 
| tkz 


If z is the polar axis, then z = r Cos® and f(r-z) has 
azimuthal symmetry. Using (2-24), the Laplacian is 
expanded in spherical polar coordinates and the 


following substitution is made. 


à - dert 9 fu) 


U= rU=Cos 09) 





с 


The substitution is carried out using the chain rule. 


2100) _ dfiw au 
ох |» | du OX 


(2-49) 


The results in the following differential equation 
for flu) 


2 o j 
E d α + (1-Uku) 4- =k | Ro=o (2-50) 


The factor k is put completely wíthin the variable 


with the transformation 


Equation (2-50) then becomes 
E J-a t (1-v) 4, T ix | Гы) = ο (2-52) 


Equation (2-52) is in the form of the confluent 
hypergeometric equation which is solved in appendix A. 
A comparison of equation (A-1) and its solution 


(A-4) with (2-52) gives 


fO) =,F Citali Iy) 
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“сг | 
Фф. = ce Feat tt Uk (r-2)) (2-55) 


A consideration of equation (A-4) shows it to be 
non-singular at the origin and therefore the solution 
де is the regular coulomb wave function. The constant 


C is determined so that asymptotically 


pem 


КЕСЕ kar 


where as üsual Q scat is the scattered outgoing 
spherical wave. In appendix A, the following asymptotic 
solution (A-19) is derived. 


‚Е Gel bl) ^- Tut aye tno Fab 


making the substitutions, equation (2-53) becomes 


‘Chabal k(r-2)) 
as ἡ 
Ф.— Ps Co) е 


cm e o ¿Cera la (e (uz y) 
Ur) a 


the logarithm is a slowly varying function with respect 


(2-54) 


to its argument. Hence at large distances from the 


scatterer, the first of these two terms is essentially 
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an outgoing plane wave and the second an outgoing 
spherical wave. The coefficient of the plane wave term 
was defined as being unity and thus the value of C 


is determined, 
-TOR f2 у 
(ze C e οὐ) (2-55) 


Equation (2-47) is also separable in spherical 
coordinates, thus leading to a decomposition of фо 


into partial waves. 
Let De = An & (eo) 


When this is substituted in (2-53) two equations for 
F and @ are obtained. 


- S. Se = + L) D =O (2-56) 


qR + (É- к. AEG = = (2-57) 


where A CL 1) is the separation constant. Equation 
(2-56) 1s Legendre!s equation which has Legendro 


polynomials as its solutions. 
№ = q (osS) 
In equation (2-51) let 


p = m E , (2-58) 


t 


With this transformation, (2-57) becomes 





2 
324 -(i fi si 52]: = © (2-59) 
Let E = al (Z) Ú (2) (2-60) 


After some manipulation, the equation for Vo becomes 


Pr rn -- (6εἰ--κλνε = © ο 


A comparison of (2-61) to (A-1) shows Va to be a 
solution of the confluent hypergeometric equation 


where, 
V) n (Rejoice | 24a le) (2-62) 
The radial coulomb wave function becomes 


LE | 
F (kr) = A € Tent! Fr(@el+cx |aera}-ackr) (2-85) 


where the Ag are determined by the same asymptotic 
considerations as before, namely that фе behaves as an 
outgoing plane wave at large r. The asymptotic 
expression for Fy is calculated in the same manner 

as for equation (2-53). By substitution of the 

proper parameters into (A-19), the following result 

is obtained. 





N 


Sl 


| І kr - 2 E 
A, €? Pat) ¿saca 


2* (2¢) PChtl—lot) = 
Al (65-92 - Le) ae! 
= а 


рса) | 


The general theory of complex variables gives 
I du 
PA+ Ei) = κ | € < 


where OF = arg( [^ (41 +ca)) | (2-65) 


with this substitution, expression (2-64) simplifies to 


E Ao à (ιο) | NET S 
с ane ES A Un e. _ бп = 


A consideration of the plane wave expansion 
in (2-17) shows the corresponding radial wave function 
to be jg(kr)., with its asymptotic form given 
in (2-19). The constants Ag are chosen so that these 


two functions approach each other as p —— oo 


A, = 24 πο ος T" 


D (2.442) 


The total solution for 0, is then the sum of all the 
partial waves. This is shown by using expressions 


(2-53) and (2-55) for Oc &nd the contour integral 





representation of the confluent hypergeometric `. 


function given in (A-7) 


= св | Е "e 
φ. = qa = & е rel, ale 


¡EE 


геу» 


a rt C Cea (c af, =) | _ 
т qu k AC (2-68) ` 


using equation (2-17), the exponential ө202(1-);; 


expanded into a series of Legendre polynomials 


where Z =r Cos8 


Do 
= ν Сб беч) e et) q (Cos 9) (2-69) 
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Where 
| = 


Е - (x8 ECT E 


e, κ. «θέ ζω πο 


a consideration of equations (2-30), (2-33), (2-57), 
(2-63), and (2-67) with equal to zero in the last 


three, gives 


X fe (х) = 247066) 00 Γ(4ε [24:2 |~2x) (2-71) 
/ PQL+) | 


using the series form of the confluent hypergeometric 


function (A-4), equation (2-71) becomes 





As 


a Ll AS a? _ 
po 28 X LEES = (2-72) 


Letting x =kr(l-t), (2-72) is substituted into 
(2-70) and the order of summation and оса 
is reversed. This gives | А 
= ICH ler 
Ë (kr) = ze (¿JŠ< e 
g ER 
ET) PC) 


x X» (04148) Ties feta prem AS Ж (2-75) 
$ 2o 


f (224245) 


The integral given in the parenthesis is a contour 
integral representation of the beta function as given 


in expression (A-6) 


ь |4 ui 7 dt ‚2 2). Poked) [ers ata) 
o [lets] 


with this О (2-73) becomes 


Elan = 24 а (куё APP tetas) 


POLER) 


x s M LIHAS) в (5645) Cade? 


x E Paetas) Pelti) : SI 
By equation (А-4), this is written. . 


Elkr) = 25897 (куе 


x Deute) F (lelece [ 2era хн) d 


r(22«2) ' 
(2-74) is compared with (2-63) and (2-67), WE 
E, ler) =. 5 (kr) e < (2-75) 
kr 
where & is defined in expression (2-55). Expressions 
(2-69) and (2-75) are combined and the partial wave. | 


expansion of Qc is obtained 


. = 7% - | 
Remind (2241)8 GO Ces) (2-76) 


Optical Wave Functions 


The eigen solution to equation (2-45a) corresponding 
to an outgoing wave. will now be obtained. It was shown 
in the previous section on coulomb waves, that if the 
nuclear potential is neglected, the equation is 
separable in spherical polar coordinates, and since Vy 
is a function of only the r variable, the separability 
is unaffected. If the Laplacian is expanded and the 
variables separated, the angular wave function is 


unchanged. The equation for the radial coordinate is 
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d* ποσο рй | a κα 
(9+ Ὄπ е) задо (ет 


The total solution is 


Fe a L = Cok LAE (G sx) ut 
x C S ΜΗ ^ 


7 


where, as usual, the constants Cg are determined 
from the asymptotic considerations. 

It was shown in equation (2-64) that the regular 
asymptotic form of the radial coulomb wave function 
could be built from a linear combination of the two 
asymptotic exponential solutions 


where the plus and minus signs correspond to outgoing 
and ingoing waves, respectively. As the nuclear 
potential was postulated as being ON these 
may also be used in determining the asymptotic form 
of Ale). This is seen by substituting the 
exponentials into equation (2-77) with the following 


result. 





Сл 
O 


. 
. 
. . 


т ла 22T +4) 


(22 + Olzal) € 


—> © 


it is seen that as long as VN AS at least 

as fast as re, an asymptotic. form of Fl kr) can also 

be formed from these exponentials with error certainly 
no greater than that inherent in the coulomb solutions. 
The phenomenological nuclear potential postulated in ae 
earlier part of this chapter meets this requirement. 
From consideration of the corresponding coulomb solution, 
the regular asymptotic solution to the radial optical 


wave function is 
fibi Sin Lkr- dlh akr -tr + 07 + Ç) (2-79) 


where όρ 15 an additional phase shift which the wave 
would experience in passing through the region of tne 
nuclear potential. The constants Cyare now determined 
from the following two considerations: (1) The nuclear 
potential, being short ranged cannot affect the incoming 
asymptotic solution and, (2) When Vy is set equal to ` 
Zero XA must reduce to the coulomb wave. As it was 
shown that the solution could be written as a linear 


combination of the ingoing and outgoing coulomb waves; 
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the first requirement becomes 
Ap Sin lke -alh akr -Lrt + ба) 
-¿(ertln2ke + б) 
e > (2-80) 
z tUer—olu2lr=2T a ox) 
+ be е > 
The solution of this equation yields 
{9% 
Ао = - 2 е 
5 2ἱ Әр | 
2 =-e l (2-81) 
u 
Са = Са < i (2-82) 
In order for (2-78) to reduce to equation (2-76) 
when Vy and thus ду are zero, 
l ¿Ls 269 
Co š = ђе (2-85) 


and the total wave function becomes 
Q) | m. Coro 
x = LY Que ZG) Plaso) (2-84) 


mere > Cm =Ta+4 éa (2-85) 





Using the Legendre addition theorem to write this 
in terms of the angular coordinates of the vectors 


X and г, à (2-84) becomes 


χ΄ἴκη - ἐπ σσ. А ¿e 


Lio mal 
x Zier) len Yes 


where the primed variables will be taken as the 


(2-86) 


coordinate corresponding to k. 
: δν 


Inelastic Integral Equation 


The following development parallels that.used 
in the a, of the general integral equation, the main 
difference being that distorted waves are used instead 
of plane waves. The inelastic component of the 
scattering has already been isolated in equation (2-43) 


and the corresponding integral equation for is 
V = —- ¿Mp He. E): lay dk: (2-87) 
EO Xo 


Kir), the Green's function is a solution to the 


following equation. Note that the f subscripts have 


been dropped as before. 


w’- 24 + е) ЖО? =-4 Уч, (2-98) 


P 





It was seen in equations (2-21) and (2-22) 
that the Dirac delta function could be expressed 


as a series of spherical harmonics as follows 


| Шм i 23 у \ 
Sie!) = L, Str) > > Yeh eea 
L: ma- l 
Equation (2-88) differs from (2-16) in the previous 
treatment only in the addition of a potential term 
which contains no angular dependence. The series 
expansion of the Green's function defined by (2-88) 


will then have the same angular dependence as (2-23). 


and can be written 


HEE) = $24 07) үө VEs (2-89) 


LO pe 


When this is inserted into (2-88) using the Laplacian 
given in (2-24), the radial Green's function is seen 
to be a solution of the equation 


2 
| < 00) EU) LU) 4 = 1 


E) (2-90) 


The reader is referred to equation (2-25) and its 
subsequent discussion, where more detail is given 


to a similar calculation. 


Кессе olor’) = r 4 (cr) (2-91) 





equation (2-90) becomes 
E l 
(ric cielo d ven 


Ihe radial Green's function is obtained from two 
independent solutions to the homogeneous part of this 
equation, one regular at the origin, and the other 
corresponding to a purely outgoing scattered wave 
satisfying the boundary condition at infinity. The 
regular solution (is obtained from equation (2-77) 
with its asymptotic form given in (2-79). 

Let Ji (ti) be the solution to (2-77) and thus to the 
homogeneous part of (2-92) with outgoing boundary 
conditions. It was shown in equations (2-80) and (2-81) 
that the asymptotic regular solutions could be written 


in terms of the two irregular exponential solutions 


gm Ln A kr AB +0) 


where the plus and minus signs correspond, respectively, 
to outgoing and incoming waves. 


Thus, asymptotically, 377 (k,1) becomes 


-a (νι ler-Lz +07, ) 
(к) — a “= di (2-93) 





The Green's function is written 


| & X) rer 
SG г) = / (2-94) 
| E Pyle) rrr! 


where, analogous to equation (2-3la), the discontinuity 


in the derivative of 9 (rr!) evaluated at r! is written 


c d Y, lkr) d | | 

€, S. Pe | e, S: Z“ F) q ( ) 
rar rest! 

The constants Є у and Go are determined from the 


Wronsklan relation 
W(Z 4) =-A (2-96) 


where AY is a constant. The Wronskian is evaluated 
using the asymptotic forms of A and 2% and found 


^ 


to be 


E d — αντ Я 
Z = Hh 2 d = -k (2-97) 
A comparison of equation (2-97) with (2-95) yields 
the values of the constants. 
| | T. i 
ET kr Ha tkr) | 
(2-98) 
E a 


2 lc t^ 


— 





Using the notation of (2-39), the radial Green's 


PA 


function becomes 
T LEN 
Zoe Z Cr) А2 (26) (2-942) 


and the total Green's function is obtained from 


(2-29) ТҮ ОО УОТ 


Hor) = aie SS yz ZK (Ey) 


z Qao NC. | (2-89a) 
X Y (o) Yet 


This is then used in expression (2-87) to obtain the 


outgoing scattered wave. 


Inelastic Scattering Amplitude 


Using (2-87), an expression for Yin the asymptotic 
region is sought. For large г, ryrt, A GJ cam be 
replaced by its asymptotic form given in (2-93) and 


the Green's σος written. 

| с (Е. —ot a2 E. ri) 
Ku ) — ак = e 
pm (2-99) 


Σο EZ OY 


waere CL, is defined in expression (2-85) and Pe is - 
taken along the direction of kp». |: 
Consider the following part of (2-99) 





O 2. 
5 ITA * ἐστ, 
г I Fi „е | (2-100) 


Upon comparison with (2-86), these are recognized 
as being similar. 

Expression (2-100) is identified with the өйы 
conjugate of no reversed" optical wave function, 
and is denoted by Ye). In the "time reversed" 
Scattering problem, the initial and final channels 
have exchanged roles ana ee) τς & solution to 
equation (2-45a) with incoming boundary conditions, 
where Vy has been replaced by its complex conjugate 
to simulate flux re-inserted into the beam. The 
CUR is perhaps best generated using the previous 


calculations for У к,г) and the Wigner Relation]? 


Hee 
| X? (kt) = ү“ -kr) (2-101) 


Using — (2-86), 


Mien = У) Мен 
x Y, (8 le, ὁ) 


The parity of the spherical harmonic is 


Yan) «αρ Y) 


(2-102) 


14 





Hence, equation (2-102) becomes 


x —) 
бр = ВОС, es 


Which agrees with (2-100) 
The Green's function, (2-99), is then re-written 


incorporating the "time reversed! wave function. 


An) 


E 
Kar & Len) (299a) 


Using equations (2-87) and (2-99a), the asymptotic 
expression for Y. is obtained. 


AN ы 


\ ) YAp (-/ 3 : 
5 M Le -= er X μέρος, (2-104) 


The first term is the elastic scattering or optical 
wave Ron which is a solution to the homogeneous 
equation (2-45) and is necessary should the initial 
and final channels coincide. This would be the case 
of elastic scattering. The inelastic scattering 
amplitude is the coefficient of the outgoing spherical 


Wave o 





И um 


With the inelastic component of the scattering 





amplitude separated, the total scattering amplitude 
can be written, not only as in equation (2-40), 

but also as a sum of its elastic and inelastic 
components. Using equations (2-42) and (2105), 


the total scattering amplitude is written, 
G Ф) 
πο = – ааг еи |197 
+ УУУУ] 


Matrix Element Formalism 


(2-106) 


Customarily, the scattering amplitude is not 
used in the literature, but rather the matrix element 
giving the probability of transition from an initial 
to some final scattering state. A consideration of the 
expressions for the scattering amplitude shows the 


matrix element to be given by 


сал 2 3 
Тег 2 +664) E^ 


where the subscripts i and f refer to the initial 
and final scattering states, respectively. 
In this formalism, equation (2-106) becomes 


the Gell-Mann-Goldberger relation . 


Ta = (Ve - {| ¿Y > UE LV! ^p» (2-108) 


which, while derived here for the specific case 
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of the optical wave functions, is quite general and is 
usually obtained from a unitary transformation of the 
general scattering amplitude matrix element, 


Te = < AF η U+V | YE (2-109) 


where Ц. and V. are any two potentials that have been 
isolated. 
The angular distribution or cross section is 


written in this formalism using equation (2-7) 


toa = Cages (D | Teal (eo 


AV 





The summation sign indicates a twofold operation; 

(1) due to the inability of the detector to distinguish 
between spin orientations, the projections of the 

final state spins aro summed, and (2) the projections 
of the initial state spins are averaged as when time 
and thus the direction of particle motion is reversed, 
the direction of any spin must also be reversed to 
preserve time reversal invariance, 1° It is seen that 
if these projections are averaged, producing ‘an ."average 
spin" initial state wave function, the problem is 
averted. A suitable modification would have to be 
made, however, if polarization experiments are to be 


considered. 
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It now Only remains to specify a suitable potential 
and scattering wave function and, using (2-110), the 
cross section can be obtained. This problem is 


considered in the following chapter. 
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CHAPTER III 
EVALUATION OF THE DIFFERENTIAL CROSS-SECTION 
FOR COLLECTIVE EXCITATIONS 


Introduction 


An exact calculation of the matrix element Tes 
is at best, a difficult problem. Transforming 
Schrodinger's equation to its integral equation 
counterpart has not given its solution, but has only 
provided a convenient way to obtain the scattering 
cross section once the wave function has been found. 
The specific interaction potential is, in general, 
not known and thus an exact solution for Y is not 
easily found, For example, if the angular momentum 
transferred in the reaction is non-zero, the interaction 
potential will not, in general, be a scalar function 
of position and normal methods for solution of the 
Schrodinger equation fail. Physically meaningful 
results, however, may be obtained with appropriate 
approximations for the wave functions and the 
interaction potential. It is the purpose of this 
chapter to discuss the techniques involved, and to 
show their applicability by fitting experimental 
angular distributions. 

In each of the expressions for Te,, four 





оК 


‘integration variables are employed. They are 

Lys Les Ts and Ep “here each is a vector quantity. 
Undoubtedly, both the initial and final position 
variables, and the initial and final internal variables 
are related, but the exact nature of this relation 

is unknown. It would be convenient to have but one 
position variable and one internal variable when 
performing these integrations, hence for want of this 


exact relation let 


їт = Ге = FE 


as 


Sr = = ΞΕ 


This is called the zero-range approximation. 


(3-1) 


The angular distribution is now evaluated using 


the general scattering matrix element. 
LE sv] ^b 
oe — U > (5-2) 


It will bo assumed that the wave function can be 


at least partially factored and takes on the form 


= ме) кю (5-5) 


where the wave function v,(&) is the eigen function 
| Hs) 
for the initial nuclear stato, and Tr, Z) is & 


solution to the equation 
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GE vu ПЕ V-E) O (3-4) 


iv should be noted that the Hamiltonian describing 

the internal motions is not Bu E Hence the = 

| LE në can be treated as constants and the function 
Bu? depends on them only in a parametric sens 

Te; is then obtained by a two step calculation; 

(1) equation (3-3), a function of only the r variable, 
is solved and the resulting solution: used to compute 


the transition amplitude for a fixed value of ©, 


mere ТЬ.) = 6 изи 949» (5-5) 


and (2) this result is multiplied by vy(%) and the 
projection of their product on the final state f is 
found. 

Physically, the above corresponds to the assumption 
that the initial nuclear state described by the function 
v,(&) changes so slowly with respect to the transit 
time ОГ the projectile, that the scattering wave function 
is completely determined by the value of vy(&) at the 
instant of the reaction. This is analogous to "stop 
motion" photography, where a moving object has its 
motion frozen on film through the use of a short duration 





flash of light. In the time considered, the E 
. variables are almost stationary and thus it is seen 
that tho Hamiltonian, HE), can be neglected in.the 
integral equation leading to (3-2) and in (3-3). 
Tnese assumptions are most valid when high energy 
projectiles are used and when the nuclear excited 
states are collective states. Since the excitation 
energy must be shared with the fairly large mass 
of the nucleus, the collective variables exhibit a 
slow change. | 

The adiabatic theory (Diffraction Scattering) 
will now be developed and, using the computed angular 


distributions, will be compared with experimental 


alpha particle scattering data. 


Diffraction Scattering 


The Fraunhofer Approximation 


The calculation is most conveniently made in the 


Scattering amplitude formalism; hence (3-5) becomes 
р | Mm 3 er 1) 
ењ е ин)? (ә) 


Equation (3-3) is multiplied through by 2m/h“ and the 
resulting expression is substituted into (3-6) thus 
transforming away the dependence on the interaction 


potential. 





o4 


Сок EZ) en 


N 
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hn E 
A^ 


} 





wnere I I 


This is put into a more convenient form by recognizing 


-ik 


that tne exponential e ar is a solution to the 


wave equation. 


If each side is multiplica by 2£ 8 and the rignt and 
left sides are, respectively, added to and subtracted 
from the integrand of (3-7), the result will leave the 
integral unchanged, but allow it to be written in 


the form . 


Farga = ge [dhe ER Qu 


ES 
(3-8) 
в ГС ερ. ARI КЕ 9 
A^ С) (те 7) 7 | 
Making the adiabatic approximation, k; κο, the first 


of these integrals becomes vanishingly small and will 


be neglected in the remainder of this development. 


Ç) 
Asymptotically Z can be written 


кг. Г | 
Le = e M AR (3-9) 
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wnere Ce: is that component of the wave due to 
‘the scattering. This is inserted into the second 


integral of (3-8) with the result 
‹ оа. 3 006-4) 6 
(как) = - 2-05-86) dre 
(3-10) 


2 ke --ἱ e = 
ο κ. αν 


The first of these integrals vanishes as in (3-8) when 
the adiabatic approximation is applied. Green's 
theorem? rolating the volume and surface integrals 


of potential functions is 


/ ФУ УФЕ = СТА tyo]. ag” (3-11) 


Vol voz 
If it is assumed that the potential U -p V is non- 
vanishing, oniy within some finite volume, then the 
second integral in (3-10) is transformed into the 


following surface integral 
L рд - ОРЕ, 
fee (Т АЕ Oy] oin 
S 


where the surface S is chosen so that the interaction 
potential vanishes outside. 

The impingent particles are assumed to be fairly 
energetic so that the high energy approximation can 


be used. In this approximation, the paths of the 
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incoming particles are assumed to follow straight 
lines parallel to the vector Κεν The surface S 
is then conveniently chosen to be a cylinder whose 


axis lies along the vector k, as shown in figure (3-1). 





FIGURE (3-1) 


With this choice of surface, TN vanishes 
everywhere except in the shadow plane in the region 
of the nuclear shadow. The entire contribution to 
the integral in (3-12) comes from this same region 


and the element of area dO€* becomes 


A 
dg- = I ЗА 
A 
where k; is a unit vector along К+ and dA an element 
of area in the shadow plane. Equation (3-12) is 


then written 


fixe =- 2 | oF (Gog, ыз (513) 


сер 
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\ 
To compute the Fraunhofer amplitude, the nucleus 
is assumed perfectly absorbing or "black" to the 
incident beam and thus ΚΘ, ο since there can be 
no particles in the shadow region. W scat - 
obtained from (3-9) by КИЕ ag o 


| t 
dium Е - e= 


This is incorporated in equation (3-11) which becomes, 


суи 
A = (14K “Gs ©) | 6 ety, (3-14) 


SHADOW 
where O, the scattering angle, is the angle between 


fu 


В. 


the vectors ki and ke. The vector difference kj kp 

is the momentum són to tho nucleus. It has 
magnitude 2k Sin ( Ф /2) if the adiabatic approximation 
| к= Ку -— is made and for small scattering angles 

has a direction very nearly perpendicular to k;. 

Thus for small O » the momentum transfer vector lies k 
in the shadow plane and equation (3-14) can be 


written 


E, OGs 


Leo E ^ rdedó (5-15) 


where the elemental area dA is written in plane 
polar coordinates, and b is the angle between the 


momentum transfer vector and re 
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The Theory of Blair and Drozdov*»% 


Tao collective variables E can only affect the 
scattering amplitude by changing the size and shape 
of the nuclear surface and thus the nuclear shadow. 
Hence the collective model proposed by Bohr and 
Mottelson is appropriate.” 

It is assumed that the nuclear surface can be 


described by the following multipole expansion 


R= > Вы SS $) (3-16) 


where the Ew are the deformation distances of the 
LM multipoles. If the incident beam direction lies 
along the z-axis, the deformation of the nuclear 
equator and thus the corresponding geometrical shadow 


can be approximated by 


p = Q, +7 8,089) | (3-17) 


LH 


This is inserted into the expression for the 
Fraunhofer amplitude given in (3-15) which becomes 


a = X Gs a) = 2 Pd 





o9 


Zquation (3-18) is written as a sum of two integrals, 
one over a disk of radius Ro and the other over tho 


deformation distances as follows; 


tee = ¿Maa J ir DE SRI E 


Zu ud 


БС (1-56) z ΠΕ ο rdrdd 


The first of these two integrals gives the 


(3-19) 


elastic scattering component and can be calculated 
by noting the following integral representation for 


the Bessel function.” 


275 


ai (2 2 
ZT die SU = | @ ssvaQdó (3-20) 
| O 
Using (3-20), it is seen that the Q integration gives 


Ф 


{= L ΠΠ J kS Onde (3-21) 


performing the r integration gives 


s tk | BICI 9, _ 
+, = 5 (1+5 $@) TIUS ο στρ (3-22) 


Using this result, the elastic angular distribution 


ís calculated from (2-17). 
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Returning to the second of the two integrals of 
(3-19), the inelastic scattering amplitude will be 
obtained where only terms linear in the deformation 
distances E yy will be retained. Integration of 
the radial variable By, parts yields. 


ры um ESAS 
ATE | 
хе Z Fimti [CS = >J (3-23) 


where A = (2k Sin. QD Cos 


e= a 


the exponential becomes, to the first order, 


E 
e £ +28 eL L| 2 ο 


When this is inserted into (3-23), the integral 
simplifies to the following. 


C2L Sii Cas d 


TE eraasju ZA. (3-25) 





To calculate the integral, the following expression 


is used? ; 


Πο ο σεξ 


QT L-M)! (М) 
(ι- μμ) (14H)! κ 


Y - 


(3-26) 


O FOR htm ODD 
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where (2n)!! = 2*4-6-*-(2n) 
and  (2n+1)!] = 1-3-5---(2n+1) 
This is substituted in (3-25), and the order of 


summation and integrations reversed. 


! lj, 
= ссе, : ? мын ea) lu) > 

ww Š ue ex ar) “un ΕΠΣΠ 
27 (3-27) 


X EIA ын, X jac C GB ὃ 0550 "Cid 


© | 
L+ M even 


A consideration of the integral shows that it can 


be written 


Tu = SERA Gs DAS y Eod Gent] > 


(L-A) (deut! 
= cit dh | Ck Ro Sn B), (5-28) 


When m is an integer, the Bessel functions of order 
m and -m are not independent and are related by 


the expression 
M 
In HUTS), 
Hence (32-27) can be written in a more convenient form 


Th = — e Ya [nu } е) je 


Zu (u)! τ; η 


tm)" (2. Yc Q. Sin) 
ΤΗ (мі 

NEM" Ἔνεκν 
The scattering amplitude is then multiplied by 





62 


v; ( £2) and the projection of this product on the 


‚ final state Vel & ) calonie can The resulting 


amplitude becomes 
fu Е (1405 9)? Gg EDS [JP 


X Tm Hl] ^ loko. SES ος 


(Lo) Cot 


The calculation of the matrix element requires 
that a nuclear model be specified. In the collective 
model, the nucleus is represented by a deformed | 
sphere of change. Excited states can be characterized 
by vibrations and rotations of this spheroid. The 
Low lying states would correspond to rotations about 
the symmety axis. 

The collective model assumes the nucleus to be 
in a deformed state characterized by a multipole 


of order L. The radius of the nuclear surface is 


then | 
ρ- Ρο(ι εάν) (5-50) 


where Al is the polar angle with respect to the 
symmetry axis and By; is the deformation parameter 
corresponding to the T multipole. By the Legendre 


addition theorem, this is expanded into the space 





fixed coordinate system. 


= Rol 1+ Bx Est] "2 


where the primes indicate the angular coordinates 

of the symmetry axis. To compute the value of the 
matrix element, the deformation operators Έτ must 
be expressed in terms of the variables used in the 
collective model. When (3-31) is compared with 


(3-16), it is seen that for a particular L state, 


Zum > (25, жа е Е (3-52) 


If it is first assumed that the nucleus is initially 
in a spin-zero state and then that the angular 
devendence of the nuclear wave functions can be 
expressed - spherical harmonics, the calculation 


of the matrix.element is simply 
<iml 100» fé E Y "V αλ. 
| 007 = (2ι-) | м μή 


The integral is the orthogonality integral for the 


spherical harmonics and has the value of unity. 


^ p, 
VIE IM? Ga ER 
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The angular distributidn of partícles 
inelastically scattered from a nucleus initially in 
a "“spin-zero" state and excited to a given multipole 


L is obtained from equations (2-7), (3-29) and (3-33). 


I 


CNAA ( 1-м)! (1+4)! 
Tor) = = ss ΕΝΝΝΝΝ ey 2 μας 


X ЈА | Heß, Sin ©) 


Ms—k, -L+2 eee 


(3-34) 


This is the basic result of the Blair-Drozdov theory. 
For the readers convenience, the elastic and inelastic 
angular distributions for the finst four excited states 
are listed below. | | 


| È ££ A Gk&,sinÉ ) 
Ute = πας, fool 


σίο»ο) - „кезй „оў? (2ER Sin) 


2 
о a 
How!) = г 2) +оо. | ek 65.6) 


| T τ > 
(>59 = LEL acosa газна) NE asit 


kel e) 


040-3) = ENTE S (rcs ST eke 5.8 OT алеу 
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The characteristic property of equation (3-34) 
is that cross sections for even values of L are 
calculated from a linear combination of the squares 
of even ordered Bessel functions, while those for 
odd values of L use Bessel functions of odd order, 
Beyond the first maximum, Bessel functions rapidly 


approach their asymptotic form.’ 


"ED @ _ to’ _. bre 
Jg ον να 2 (3-35). 


(Tk Ro Sin OL >) !⁄> 


Hence it is seen that the angular distributions for 
all even parity transitions have the same shape. 

It follows, then, that all odd parity transitions 
have the same shape. The two, however, oscillate 
exactly 180 degrees out of phase with each other. 
It is further noted that the angular distributions 
for all odd parity transitions are exactly in phase 
with the elastic cross πέσο η thus providing a 
means for parity identification. This rule, known 
as the SIEHE rule, is usually well obeyed in 


experiment. 


Analysis of Data 


1 


A consideration of equation (3-34) shows that 
there are two adjustable parameters within the context 


of the theory; Rọ» a quantity of the order of the 
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nuclear radius and/O,, the deformation parameter. 
Rg is found by fitting equation (3-22) to the elastic 
cross section. | 

Fits to ara, arte data are shown ín figures 
(3-2) and (3-3) for incident energies of 32.8 and 
41.0 Mev(Lab), respectively. R, was taken as 6.45 
Fermis in both cases. For those angles considered, 
the Blair theory predicts the oscillatory nature of the 
cross section, reproducing the maxima and minima very 
well. It does, however, unlike the experimental 
distribution, have zeroes. Іп the back nass, the 
magnitude of the theoretical prediction drops off 
more slowly than the data and fits, in general, become 
much worse. 

It should be noted that due to the large 
MU ions of the elastic cross section, these figures 
show plots of the ratio to Rutherford scattering, 


where the Rutherford cross section? is given by 


05 = = ES Coe © 


With the radius parameter Ко determined, the 
inelastic distribution may be calculated to within 
a normalization constant using expression (3-34). 


The theoretical distribution is then normalized to the 


experimental cross section and Puis obtained. 
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There is some concern, however, as to which peak 

is to be fitted since it is expected that the Blair 
formula will be inaccurate at small scattering 

angles (<20°) because of its neglect of coulomb 
effects and at large angles ( 390°) because of various 
small angle approximations. 

In the derivation, it was assumed that particles 
having impact parameters less then Ro would be absorbed, 
and those with larger impact parameters would pass 
unaffected, that is they would not be scattered. For 
charged projectiles this is an obvious oversimplifi- 
cation, as certainly two charged particles passing at 
distances of the order of the nuclear radius would 
experience a deflection. 

The apsidal distance, or distance of closest 
КО сасы, d, of a charged particle scattering from 


the coulomb field of a nucleus is given by? 


2 
CS HE (+ sc 2) 


where Z and Gv are respectively, the charge numbers 
of the target and projectile nuclei, E the center of 
mass energy of the incident projectile and Ө, the 
center of mass scattering angle. This is solved 


for ð, giving, 
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@ = | 


if the apsidal distance is taken to be the order 


of Ro» then for 41 Mev alpha particles incident 


36 


on Ar””, the scattering angle is of the order of. 


22 degrees. The corresponding result for the 32.8 Mev 
alphas is 27 degrees. The "black" nature of the 
nucleus will not allow impact parameters less then 

Ко and thus the coulomb scattering should fall off 
beyond these angles. Hence normalization of the 

cross sections to the second peak seems justified. 


10 


The level scheme for the first three low 


56 


lying states of Ar”” is shown in figure (3-4). 





FIGURE (3-4) 


Fits to the ar95( to!) ar? даба are shown in 
figures (3-5) through (3-8) for 32.8 Mev and 41.0 
Mev alpha particles. It will be noted that the 
product &Ro is more customarily quoted in the 


literature, and will be used here. The first 





τν 


excited state is strongly excited and it seen to be 

a 2°” state. Hence it is fitted using equation (3-34) 
with L set equal to 2. The results for the lower 
energy data are shown in figure (3-5) where /5Ң = «82 
Fermis and for the higher energy data in figure (3-6), 
Where {oRo 2 +71 Fermis. It is interesting to note 
how well the Blair phase rule is obeyed. A comparison 
of the La 2 data with the corresponding elastic data 
shows them to oscillate exactly 180 degrees out of | 
phase. 

The second excited state is strongly excited, 
with mes spin known to be either 2 or 3. Upon 
comparison of the experimental distributions shown 
in figures (3-7) and (3-8), with the corresponding 
elastic data, they are seen to oscillate in phase. 

If the Blair phase rule is correct, then the second 
excited state can only be a 3 state and thus its spin 
and parity are determined. The angular distributions 
are calculated using (3-32) with L set equal to 3. 
These results are also shown in figures (3-7) and 
(3-8) where „Н Ξ +67 Fermis for the lower energy 
data, and .56 Fermis for the higher energy data. 

The general appearance of the inelastic angular 
distributions is similar to their elastic counterparts. 
The о os nature is predicted and the position 


of the maxima and minima are reproduced well. 
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As in the elastic case, the inelastic distribution 
predicts zero minima and fares poorly in the back 
angles. 


e 


The Distorted Wave Born Approximation 


This and the remaining text follow the treatments 


li ana McCarthy. 


given by Austern 

If it is assumed that the average interaction 
potential U(r) of the colliding particle is primarily 
responsible for determining the optical wave functions 
X. (cr), and if Кү < Кс, then the potential V (m, E) 
may be treated as a perturbation. Consider the 


matrix element Tp, giving the inelastic component of 


the scattering. 


Te = CU ИГУ (3-36) 


The wave function ^P is then approximated by the 
outgoing optical wave function that was used in 


(2-41) to compute the elastic scattering component. 


р = моку. (әт) 


where the approximation of (3-1) has been applied. 
With this expression for 3D, equation (3-36) becomes 
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3p. а V (3-38) 


The width factor vgl 1: influences the magnitude 
of the cross section and depends on a detailed 
knowledge of both the reaction processes and the 
nuclear wave functions. The distorted wave functions 
ps (k,r), are primarily responsible for determining 
the angular distribution. | 

Ihe potential VeC Oo €), assumed er function 
of position, is expanded in multipoles of the vector r. 


This expansion becomes 


34 
УЕ) = 2 Vigne) Yee) κ (3-39) 


where it is necessary for the ensuing calculations 
that Уту behave like the spherical harmonics under 
rotations of the coordinate system and have parity 


(ec The factor Lv Vr Ivi > can then be written. 


α΄ 
<V. l V iv? =) C. I Vul Ve? M. (3-40) 
LM 


‘The wave functions vp and vy are assumed to: 
have their angular dependence described by the 


spherical harmonics. They will, in general, be a 
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function of their angular momentum J& and its 
projection Ma. A fundamental property of the matrix 
elements of these spherical harmonics and functions 
which transform in & like manner under'a coordinate 
rotation, is the Wigner-Echart theorem, 1? 

Applied to the above "width" factor, it is 


SARUM TAA s QLMEH JUR) VID (er 


Js and Je are the initial and fínal nuclear spins 

with М, апа М. their respective projections. The 
factor σταρ l Vy l КО» is the "reduced width" and 
(J,L M.M | JpMp) a Clebsch-Gordan coefficient. The 
reduced width is just the width factor expressed in 

a diagonalized representation. Using (3-41), equation 


(3-40) becomes 


«у. ММУ =) (SL MLM |2,M,)<2 MM) Y. (3-42) 
LM | 


The interpretation being that the LA multipole їй 


equation (3-39) transfers an amount of angular momentun 
L to the nucleus. The Clebsch-Gordan coefficients 


are so designed as to limit L to the range 


| x= Jy | €L<{Jr+Je] 
and thus insure conservation of angular momentum. 


The "reduced width" factor is a function of only the 
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r variable and it becomes convenient to write it 
as the product c. a "strength factor" and a "form 


i actor" which gi" s its radial dependence. 
№. МУ = А, Б el 


The strength factor A, would give the magnitude of 
the reaction and requires a specific knowledge of both 
the nuclear wave functions and the nature of the '. к. 
reaction. The form factor, however, need not be 
known exactly, and meaningful results are obtained 
from knowing its general trend as will be seen later. 
Upon substitution of expressions (3-42) and 
(3-43) into (3-38) the matrix element becomes 


Бел FA (AE М.Н E Mp) Ra u=: i | 
LA 
where | (3-44) 


NE = [45 XP ТҮ уу” 


The cross section is then obtained from expression 


(2-110), 


el: Oy г 2. | ən AL 


P 


(3-45) 

QJ; L MH [Je He) σι 

lt will be Са 15 that the summation sign referred 
to a sum over Mp and an average of My. For an 


initial state having spin Tis there are (2513 1) values 
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of M;, hence (3-45) becomes 


Lar σπα a 2) +1 2222 кё 


Mp AM рм! (3-46) 
X GEL MEM MOLA 10, Hr) | 2 


Noting the following Clebsch-Gordan symmetry 
relationship ; 


кй 
(м,н |.» Μα) z Cl) ИСЭ ATEO (3-47) 


the cross section can be written 


pu = с Sp (23540) SEE 2 2 
GR S 3 (23:4) 2 2 IA, ы 

ais ΗΝ à 
x LIRA Е И Je J pH Mel LM Ozde-Me Em) 


The Clebsch-Gordan ЕИ exhibit the following 
ORTMOGOM ALL 15 





(3-48) 


NE e 
> I 

L CizJe-MxMel LA Y (eJg -MpM, | bm) = б "C 
Fre 


Hence the cross section expression simplifies and 


is written 
mue (Ice Je +! x PS vr 2 | 
TUE ta : 
στ απ τε еч). буз 7_| d (3-50) 


If the spin of the initial state is zero, then the 


— 
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final state spin is just the angular momentum 
transferred in the reaction. It should be noticed 
that different values of the momentum transfer L 
do not interfere in determining the reaction. 

The problem of evaluating the cross section 
has now reduced to that of calculating the matrix 
element Č rx- llora e ¿rita “Orms Of the optical wave 
functions are required. These may be approximations >>?” 
or with complete rigor, the spherical harmonic 
expansions derived in chapter II may be used. The 


following sections discuss the nature of possible 


approximations. 


E Plane Waves 


Wnile it may seem paneer to include a 
section with this title in a discussion of the 
distorted wave Born approximation, the use of plane 
waves is instructive as it gives insight into the 
nature of the distortion necessary to reproduce the 
form of the experimental angular бузы К σας. 

It is convenient to make the identification 


Sir). JE 
ptr) a = x. 2 (3-51) 


with the matrix element Ç Ім given in (3-45) 


then becoming 


£M 





(Ste) M 
MN = [pore 7 |, ed d o^ (3-52) 


The quantity (г) o19 Un). a5 nenas ku i 

and may be thought of as describing the properties 

of a probe which measures the evolution of the nuclear 
state during the process of inelastic scattering. 

To clarify this idea, consider the case where the 
Optical wave functions have been replaced by plane 


Waves. 
де?” - SEL (5-55) 

where K = k; = Kr and Po is a constant. Here the | 
probe measures the difference in the wave number 
vectors for the initial and final states and the 
matrix element Yyy is proportional to the probability 
that the initial and final states differ by momentum 
КА. Since the inelastic scattering cross section is 
obtained from this matrix element, the use of plane 
waves implies. that the entire momentum transferred 
goes into the inelastic component of the scattering. 

If / (x) is anisotropic in space, the probe tries 


to make a simultaneous measurement of its position 


the resolution of the momentum transfer measurement 
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‘and the momentum transfer. By the uncertainty principle, 
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must then be reduced, In putiente > theory of nuclear 
reactions, it is assumed that the reaction is confined 
to the surface. It would then seem that the probe 
could obtain no information whatsoever about the 
radial component of the momentum transfer. On the 
other hand, since no angular localization is assumed 
on the nuclear surface, the resolution of the angular 
momentum transfer should be a maximum. This is seen 
in the following calculation. | 

Setting the constant Po equal to unity, tne 


matrix element C ry becomes 


Ge = е“ —— DA Ót-- e) ed (3-54) 


where F(r) becomes a Dirac delta function limiting 
the E to a spherical shell of radius Hg. 
Using eque (2-18), the plane wave is expanded 


in a series of spherical harmonics and (3-54) becomes 


/ ο, | | 
1, 2 6 uten) Lr] ο (3-55) 
Li‘ 


The integral in (3-55) is the orthogonality integral 


for the spherical harmonics where 


P 
VIE! = δι,’ Sum‘ (3-56) 





e— = == Rau - E u m 
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The Kronecker delta functions are zero unless L' and 
Wit equal L and M respectively, hence only the IM 
term remains. | 


Equation (3-55) becomes 


: | 
Tu = Z (kl) У ее) (3-57) 


where the primes indicate the angular coordinates 

of the symmetry axis. As it will be assumed that the 
nucleus is initially in a spin zero state, the cross 

section for a given final state L is given by (3-50) 


to be proportional to 


- 


€ 


L 12 
δ; ος | 4 Аке, |^ 2 1.80) [ cos) (3-58) 


CN 


The Legendre addition theorem is used to convert the 
sum of spherical harmonics into & Legendre polynomial 
of argument unity which has magnitude unity for all 
values of L, Equation (3-58) veio 


(3-59) 


O] a Ai (kee) 
Beyond their first maximum, these spherical 
Bessel functions rapidly approach their asymptotic 
form given in equation (2-19), and thus in this region, 


the cross section can be written, 
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0 —- Sın (KRo- N (5-60) 
e 

Hence this behaves like the Blair formula. It 
appears then, that the sharpness of the maxima. and the 
fact that the minima actually reach zero can be 
associated with the maximum resolution of the angular 
momentum transfer. If this interpretation is correct, 
it should be possible to spoil the resolution by 
causing /J(r) to be non-uniform over the nuclear 
surface, 

We proceed by localizing the reaction to part 
of the nuclear surface. As alpha particles are 
strongly absorbed, those incident on the nuclear 
surface will have their wave funcions considerably 
reduced on the shadow side of the nucleus: as compared 
with the illuminated side. The same effect would be 
true of the scattered particles. The overlap of ree 
wave functions would be greatest on the side of the 
nucleus opposite to the direction of momentum transfer 
Piet а shown in figure (3-9). It is assumed that 
that part of the surface axially symmetric about the 
direction of momentum transfer and cut off by a cone 
of semiangle œŒ contributes uniformly to the matrix 


element, and the remaining part of the surface 
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FIGURE (3-9) 


, contributes nothing at all. The matrix element l im 
is calculated using the expression in (3-54) where 
the plane wave has been expanded in a series of 
Legendre polynomials. After the radial integration, 
(3-54) becomes 


и“ 2. сё (рм) Z Ga) 


2/8 el 


ш Л (64) Pitos6)Sinedadp \°™ 


where the axis of symmetry is taken along the direction 
of momentum transfer. The Legendre polynomials and 


spherical harmonics are related through the expression 


РС»). £ Yee) (3-62) 


and as the azimuthal integration has been unaffected, 
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it is seen by (3-56) that the integral is non-zero 
only when M equals zero. Using this result, the 


matrix element Cry becomes 


| οἱ 
iz / a^. 
Un ος 6 Cut) ote | god G ado (3-63) 
ТАЕ 
| | / 


Ihe integrand is expanded in a series of 


Legendre polynomials. 


Yeo) Ya (5$) dh ο, 7 71 (ces) 


- Z, ë Yeso) 


The expansion coefficient is obtained by multiplying 


(3-64) 


both sides of (3-64) by Y, (6,0 42. ana integrating 
over 4U7 radians. The orthogonality of the spherical 
harmonics determines the ан "LL x on the 


right oe giving the expression 


The following relation involving Clebsch-Gordan 


.coefficionts is noted, 1? 


AL, е 7 Ufo 
Nauen 226527 


x Ch 15 «νά. αν ο» 


(3-66) 





Using this relation, equation (3-65) becomes 


(2241) (2241) ] a 2 
C NE mab (LLoo llo) (3-67) 


Ihe expansion in (3-64) is then obtained. 


le] {ο | | 
Ihe e = rar J Í (1100140) 2 (cos 0) (3-58) 
Ё. ғо 


The matrix element тз becomes 


(dtm of Z οἱ (cs) u (66, (4400 Ho)* 
ο x 
X / 4 (&:6) Sedo 


(3-69) 


Hence for a given angular momentum transfer L, the 
matrix element is proportional to a sum of spherical 
Bessel functions of different orders. In general 
this will be complex and it is therefore highly | 
improbable that Zr will go to zero for any en, 
value of K., The angular rossini ОО синд а with 
zero minima has been spoiled. If Oz 7C, the entire 
surface then contributes uniformly. The integral 

in (3-69) becomes the orthogonality integral for the 
Legendre polynomials, giving zero unless 2=0. 

With «0, the Clebsch-Gordan coefficient will be 
non-zero only when L' = L, hence it is seen that the 


result reduces to that in equation (3-59). 
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An interesting interpretation of (3-69) is 
that for total angular momentum transfer L, L! 
represents that going into the inelastic component 


of the scattering as it results from the ο. 


partial 
wave in the expansion of a plane wave, and & that which 
has gone into other components, such as elastic 


scattering or absorption from the beam. 


Approximate Optical Wave Functions 


Based on a study of optical model wave functions, “97 "+ 


reasonable approximations“ for the incoming and 


outgoing wave functions are found to be 
e л д 
Же = ΑΘ 
сю AA 
# heb tbr) (5-70) 
XO x AG) Ó ^^ g 


where y is referred to as the anisotropy parameter 

and is a measure of the non uniformity of the nuclear 
surface. When these expressions are used in equations 
(3-51) and (3-52), the matrix element Cry is seen 

to We | | 


(5-71) 


Ä EL 
σι, = E [09 £0 [AU ο 


Where as usual, K = ky ~- kp. Let, 





£ g (4.3) (3~72) 


and expand the resulting plane wave in a series of 
spherical harmonics. The angle integrations were 
computed in equations (3-55) and (3-56) and they 


result in the following expression for (a 
> 


x 
Go = 2:29 (б, afia 6) ο (373) 
д 


The cross section is calculated using expression 
(3-50), where the spin of the initial state is assumed 


to be zero. 


PF. 
cxc2 ped YR - ena) 


x | / Eco Mol zer car 


Using the Legendre summation formula, 


L. 4. | 
An? L9 ln By zt) (Cost) (3-75) 


< 
where [|' is the angle between Z ana Z ° 
СозО) сап be written 


F, = 2 


35% 


(3-76) 
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If the adiabatic approximation k; © kp» =k, is 
employed, 5 can be written 


" | TENA 2-77 
3 == 2 (K+ of) Sin S & ( ) 


where @ is the amgle between Ks, and Kp anaes is а 

unit vector along the direction of momentum transfer. 
(5-75) is inserted in (3-74) and it can be shown that 
Соз( П.) = 1, Hence Pp, (Cos Г) = 1 and the cross 
section for a given value of L is proportional to 


« 
(3-78) 





a «| f ce) Vaca, Ge Pee 


lt now remains to make some approximation for 
the radial integral. Since alpha particles are 
bns οι by the nucleus, | alr) | š will 
fall off for r less than the nuclear radius. On the 
other e F(r) being associated with the nuclear 
wave function, will fall off for r greater than the 
nuclear radius. Assuming that F(r) is slowly varying 
over the mean path of the alpha particle, the radial 
integration should essentially depend on the radius R. 
and the thickness of the surface. A simple empirical 
choice which involves only these two parameters is 


. 2 
2 — (rn ё.) 
ACI 2 (3-79) 
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The final form for the differential cross section 


is therefore 
= ame 
O J e A, (gr) rd (3-80) 
© 
With the adiabatic approximation for Z given in 


(3-77), (3-80) can be written in a simpler form. 


Z | 2 
OL £L О e Al (кР сі) Sih [rer (3-81) 
Q 





The corresponding surface interaction result is 
ә " * o A 
0, « | (> [ al is Re -cr) ma (2-92) 


In the discussion resulting in equation (3-69) it 
was found that by introducing an angular anisotrophy 
into the mm the resolution for a measurement 
of the angular momentum transferred was at least 
partially spoiled. 

It was seen that the matrix element, Cm: was 
no longer proportional to a single spherical Bessel 
function of order L, but consisted of a sum of е) 
The result-in (3-82) is but a single spherical Bessel 
function. However its argument is complex. That this 
actually corresponds to (3-69) is shown by expanding 
(3-82) in a series of j;' (2kRoSin@). 
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Consider the plane wave expansion, where the 


magnitude of KR, = x +> 
č (x+y) Cose p 
= "2t (2ne) 4, (xey) Ê > (ese) (3-83) 
This can also be written 


CX Cv 4O wy Cos A 
Q ЕС 5 (244) (4) 
Б `| 


αν Ú| (y) Р, (Со; ©) P (Coso) 


(3-84) 


The two Sondus are equated and. both sides multiplied 
by P;(Cose) Sinodo. When the result is integrated 
from oto 7C, orthogonality eliminates all but the 


p^ term in (3-83) giving the following expression 


" Key) = уема (24%/)(24+/) 144. (4) 
vt (3-85) 
х2 0000) Фо) 0з ө) nodo. 
The integral is written in terms of spherical harmonics 


and both sides are multiplied by db and integrated 
from O to 21 T 


4 Gu) = 2 2f" M ~) 24 M .. 4) 


(3-86) 
M > / x 
SA, n 
X N Tan) (2041) =D A Jo, d. 


The integral is given in (3-66), and with the proper 
substitutions (3-86) becomes 
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š ——À | | 
j. 6e») =? ALL e ibo Me) ful dat) (3-87) 


d 
BG 


Letting, 


X= IEE 22 @ 
у= асі 5» @ 


It is seen that. 


4, L2 (eg, CY) Sin 27 = 


¿tes ^ 1) 


| , T (3-88) 
X (LLoo/to) Je (2c y Sin£ ) £ (26е, Siu & ) 


which is the desired result. The parameter Y 
evidently has the same effect on reducing the resolution 
of an angular momentum transfer measurement as the 


introduction of an angular anisotropy into the reaction. 


Results 


Optical model fits to the ar°°(a,q)Ar™ data 
are shown, for the lower energy case in figure (3-10), 
and in (3-11) for the higher energy case. These are 
included as the approximations in (3-68) are based on 
a study of the optical wave functions.and thus the | 
ability of the optical model to predict the elastic 
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distributions should indicate the validity of these 
approximations for the data considered. The maxima 
and minima are reproduced very well, both in position 
as in the Blair formula, and in magnitude. Although 
not indicated, the optical model does well in 
predicting the behavior of the back angles. A radius 
of 5.6 Fermis was used in both cases. 

Fits to the Ar°°( qm, Atar? data obtained using 
(3-81) and the appropriate choice of L are shown in 
figures (3-12) through (3-15). А single set of the 
parameters, A and Y were used to fit these data. 
They were As 0.88 and Y/2 0.90 as was suggested by 
McCarthy and Pursey in the text of their articie.. 
Similarly a single value of the radius parameter В, 
was used. This was/found to be 7.0 Fermis. It is 
seen that the angular predictions of. the maxima and 
minima are comparable to the Blair theory». There is, 
however, improved agreement in the relative magnitudes 
of the theoretical and experimental distributions. 

The non-zero minima associated with the poorer 
resolution of the angular momentum transfer measurement 
is clearly shown. A systematic study of the parameters 
À ana Y was not made. The effect of Y, however, 

can be inferred from the previous study. Setting Y 
equal to zero removes the anisotropy and thus the 


minima will again reach zero. The interpretation 
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of the series in (3-88) was that the entire angular 
momentum transfer Lh no longer went into the inelastic 
component of the scattering but only an amount L'À. 

The rest was absorbed out of the beam or went into 

the elastic scattering component. Hence it is expected 
that Y would also affect the rate at which the cross 


section decreased with angle. 





MEME ; 


100 


DIFFERENTIAL CROSS SECTION (MB/STER) 


99 


Айх AGO”. 
Q =-1.98 MEV 


Ex = 32.8 MEV(Las) 


10 20 50 40 50 60 70 со 99 
CENTER OF MASS SCATTERING ANGLE 


FIGURE (3-12) 





DIFFERENTIAL CROSS SECTION (MB/STER) 


100 


100 


Aria Ajos 
О =-1.98 МЕМ 


Ex = 41.0 ME Vías) 


10 20 30 40 59 GO 70 CO 90 
CENTER OF MASS SCATTERING ANOLE 


FIGURE (3-13) 





100 ; 


(MB/STER) 


DIFFERENTIAL CROSS SECTION 


20 


Ач ,«)да?* 


Q=-4.17MEV 


Eu = 32.8 MEV (LAS) 


$0 40 50 co rae) CO e 
CENTER OF MASS SCATTERING ANGLE 


FIGURE (3-14) 


101 








100 


DIFFERENTIAL CROSS SECTION (MB/STER) 


IO 


А Да 
Q=-4.17 MEV 


Eu = 41.0 MEV(LAB) 


$$ 


20 50 40 50 GO 70 со со 
CENTER OF MASS SCATTERING ANGLE 


FIGURE (3-15) 


102 





103 


NOTES: CHAPTER III 


1 Walsh, John B., Electromagnetic Theory and Engineer- 
ing Applications, (Ronald Press, 1960), p. 36. 


2 Blair, J. S., "Inelastic Diffraction Scattering, " 
Physical Review, 115, 928(1959). 


К GT Drozdov Sal "Ihe Scattering of fast Neutrons 
y non-spherical Nuclei," J. Exper. Theoret. Phys. USSR, 
28, (1955), p. "734, p. 736.. 


4 Bohr, A., Mottelson, Ben R., "Collective and Indivi- 
dual-Particle Aspects of Nuclear Structure," Det Kel. 
Danske Videnok. Selskab, Mat.-fsy. Medd., 27(1955). 


o Morse, P. M., Feshback, H., Methods of Theoretical 
Physics 2, (2 vols., McGraw-Hill, 1955), р. 1525. 


6 Adler, Bohr, Huus, Mottelson, and Winther, Reviews 
of Modern Physics 28, 432(1956), Eq. II A. 25. 


7 Margenau, Henery, Murphy, Go Moe, The Mathematics 
of Physics and Chemistry, (D. Van Nostrand Company, 1956). 


8 Dicke, R. H., Wittke, J. P., Introduction to 
Quantum Mechanics, (Addison-Wesley, 1960), p. 61. 





9 Becker, Re Ae, Introduction to Theoretical Mechanics, 
(McGraw-Hill, 1954), p. 240. - 


10 Endt, P. M., Van Der Leum, C., "Energy Levels of 
light Nuclei. III," Nuclear Physics, 54, (1962), р. 216. 


ll Austern, N., Selected Topics in Nuclear Theory, 
(Vienna, 1963). 





12 McCarthy, Ie E., Pursey, De Lo, "Simple Realistic 
Treatment of Nuclear Direct Interaction Processes," 
Physical Review, 122, 578(1961). 


13 Preston, M. Ae, Physics of the Nucleus, (Addison- 
Wesley, 1962), Eq. A-"7S. - 


14 Lbide, Ба, Α- ὅ]ο 








104 
16 Rost, E., Austern, N., "Inelastic Diffraction 
Scattering," Physical Review, 120, 1375(1960). 


17 Marr, Ge, Unpublished M.A. Thesis, Miami University, 
Oxford, Ohio, 1965, 


18 Butler, S. T., "Direct Nuclear Reactions," Physical 
Review, 106, 272(1957). 


19 Preston, Physics of the Nucleus, Eq. A-69, 


20 Eisberg, R. τως осер lo Es», Nuclear рос 
TOO 7 1 (2059 0 = 





оса, Т, Ee, Nuclear Physics, 10, 583(1959), 
574 (1959). | 


22 McCarthy, Pursey, Physical Review, 122, (1961) 
Pa 98l., 


\ 





105 


APPENDIX A 
THE CONFLUENT HYPERGEOMETRIC EQUATION 


Kummer 's equation κε hypergeometric equation 
often encountered by physicists in the solution of 
Schrodinger's equation. 

It has the form 


d'u ς _ = 
298+ (-2)94-94=0 um 


Solutions of this equation are known as confluent 
hypergeometric functions. This appendix will present 
only enough of the theory to provide an understanding 
of chapter II.: The reader is referred to references 
/71_7 апа / 2_/ for a more complete development. 


The Series Solution 


It is assumed that equation (A-1) has & solution 


OO 

п 

q= È Cn“ 
X 


This is inserted into the differential equation 


of the form 


resulting in & polynomial in z. For y to bo a solution 
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Por all values of z, the era of all powers 
of z must be identically zero. With n set equal to 
zero, the coefficient of the lowest remaining EE 
of z yields the indicial equation. 
Solving this, gives 

k =0, l-b (А-2) 
The first of these gives the desired solution which is 
regular at Эте origin. The following recurrence 
relationship for the coefficients is then obtained, 


assuming that ορ is non-zero. 


С = А _ 


(n +1) (b + hn) 


With this, all the coefficients can be generated 


(A-3) 


to within an arbitrary constant c, and setting this 


O 
equal to unity yields the following series : 


D (an) IO) 27 = = 4 
Ч = 2 πα) F (p+n) nl = ,E (al biz) A-4 


À Contour Integral Representation 


The expression (A-4) is re-written 


a (6) Parn) h-a) z α) κ΄ 
e (aibl2) Pa) 0) AE P (bn) ni 


Fco) Tatas) Boon bd d 2 (А-5) 


P G3 Co-2)) 





107 


where B(u,v) is the beta function and is defined? 


l 
(чу [' (v) wf и! 
Blu,v) = Fu) Pi. t“ (1-6) dc 
Г(ц+у) 
O 
A contour integral representation of the beta function 


can be shown to be 


В(им) = (1-ё Ре i |- ) dt (A-6) 
C5 
where the contour C, is described in figure (A-1) and 
the phase defined, so that along the segment from 
l to 2, 


li 
о 


КЕЕ (b) = arg (1-5) 





t=0 vd 


FIGURE (A-1) 


For the Medion to take on the same value after a 
complete circuit of the contour, it is seen that u + v 
must equal an integer. Equation (A-6) is inserted 
into (A-5) giving 


geile) a 


F (aliz) = DOUE E £P gue ας 
T'ta3T (bo) "m A 





G 
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Interchanging the order of summation and integration, 
where it is seen that b must be an integer to satisfy 
the requirement that (A-6) be d valued and thus 


(А-7) be a solution to (A-1), En becomes 


Tei (b-a) 5 —/ 
PÉ τι e al (y, eo^ 
Б. (0162) = Para) po ni ¿E -έ 
(A-7) 
r (5) (1- Uem шо ie ye “Ye 
rayito l'Cb-4) 


When the integrand of n is substituted back 
into the differential equation, the result can be 
expressed as a perfect differential and thus any contour ' 
which has the same value at its end points for which | 
the integral converges leads to a solution. In 
particular, two of these will be considered yielding 


irregular solutions of (A=1). 


e Ка) = МЫ) ея) оа) 


where 
rt (1- e ау! at g- дё 
υπ ο = pia e t (-t) "de wv 


Cj Co. 
Ihe contours C4 and Co are described in figure 


(A-2) 
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FIGURE (А-2) 


where for convergence of the integrals, ' 


1/2 < 449 2 + 172077 < zz 
is the argument of the end points at infinity. 
Consider the expression (A-9). By convention, 


the substitution zt = -u is made and 


| шщ a ua, „pa 
TOLEDO 2 24040) de 
Wileldiz) = rare) Г C 2 Li 


Ф. 


where the contour in the U-plane is shown in figure 


(А-5) 


λ 
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FIGURE (A-3) 
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Here, arg (u) = О а1опд the bottom portion of the curve 
along the real axis as was defined in figure (A-1). 
With this contour, 


RK arg z < πα 


Letting first t=t +1 and then zt =-u, an expression 
for Wolalb]z) is obtained 


52004 yola + α-ί 

V. (Ыг) = YETI) SR me _y\,, b-a-! 

š | a.) Ct -a) e ( £ u А“ (4-11) 
LU | 


The contour C. is shown in figure (A-4) 


qrg () => 


(2.20 
FIGURE (А-4) 


where arg(u) = 0 on the top portion of the curve along 
the real axis, again to be consistent with the phase 
convention in figure (A-1). To obtain convergence 


Asymptotic Expansion 


Asymptotic series expressions will be obtained 
for the irregular solutions Wa end Wo and from these 
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and (A-8), an asymptotic expression for ‚F7lalblz) 
will be found, 
Consider expression (A-1210) and its associated 


contour. The binomial expansion is used to obtain 


— n 
(i+ мун. n(b-e) B 
z ü U(e-a-«y vl 


hao 


This is a series of negative powers of z and is thus 
an expansion for large z. This is inserted into 
(A-10) and the order of summation and integration 


reversed. 


[Ыгу = ри) E ПАМ | 
W, (alba an 5 (1-6) 


Aso 
= E t 

mo Г pas nt 
C4 


The remaining integral and its associated contour 
C, given in figure (A-3) can be shown to be a 


contour integral representation of the gamma function. 


de“ ny S (A Plan! 
C4 


Equation (A-11) becomes 


-— nu er 5- Pla +n)p (b-9) = (A-12) 


τον P (6 -a.-n) du ni 


This is put in a form where (<2) is the variable 





TI 


of tne series with the following identity which can be 
verified by induction. | 


DUC) pp ya Dan) В 
Cann) σι) "πο. cin 


where n is an integer. Equation (A-12) becomes 


W, (alolz) = 





E 6.) Plants) Ca)" rr 


c nla) Γία-ὐ-εἰ) Ὡ {. 


In expression (A-11), the binomial theorem is used 


to obtain 
ë? Пса) шу 
(i- uf» 2-8 Tla-n\n! 


Equation (A-11) becomes 


| і, - | 
| ; KOLIS —]) α-ξ 
We bla?) Pa) ce) © 
(A-15) 
Pb) Gon, i ha 
x2 A h! 


As it was required that b be an integer, this can 


be written 
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- 
(alblz Tes р Nla) (-2) 
ο» - 2 , P(C&-h) PQb-s)n! 


| | (А-16) 
X (g^ "5-9 у" ра а and Ne 


Ce 


Ihe remaining integral and its related contour can be 
Shown to be a contour integral representation of the 


gamna function 


Peru du ao, 1) Toman) ` 
Ὃν 


Equation (A-16) becomes 


6) 2-55; — D) P (bo. «n) C7. (-2)” τ _ 
Walbla) < e? E Бл. (4-17) 


. Letting a —--a in (A-13), equation (A-17) can be 


transformed in a manner similar to (A-14) 


Di = Us 
Wo ell ч. To. ds ‚> TS cipal : Jr (1-18) 


For very large z, only the first term of the 
expansion is retained, thus equation ao gives for 


the regular asymptotic solution 





114 


| "m Lu > Br Be 2 a-b (A-19) 


where 
9а < тз < 
in the development used in chapter II, z takes on 


pure imaginary values, `i% 


arg(z) = + Te 
In order to use expression (A-19) A limiting procedure 
will have to be conceived where arg(z) is allowed to 
approach 7/2 or -7/2 in the direction of positively 
or с increasing angles, Ξ-- Lt will 


be assumed that this is done in all expressions therein. 
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